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The PLATO satellite mission project is a next generation ESA Cosmic Vision satellite project dedicated to the detection 
of exo-planets and to asteroseismology of their host- stars using ultra-high precision photometry. The main goal of the 
PLATO mission is to provide a full statistical analysis of exo-planetary systems around stars that are bright and close 
enough for detailed follow-up studies. Many aspects concerning the design trade-off of a space-based instrument and its 
performance can best be tackled through realistic simulations of the expected observations. The complex interplay of 
various noise sources in the course of the observations made such simulations an indispensable part of the assessment 
study of the PLATO Payload Consortium. We created an end-to-end CCD simulation software-tool, dubbed PLATOSim, 
which simulates photometric time-series of CCD images by including realistic models of the CCD and its electronics, the 
telescope optics, the stellar field, the pointing uncertainty of the satellite (or Attitude Control System [ACS] jitter), and 
all important natural noise sources. The main questions that were addressed with this simulator were the noise properties 
of different photometric algorithms, the selection of the optical design, the allowable jitter amplitude, and the expected 
noise budget of light-curves as a function of the stellar magnitude for different parameter conditions. The results of our 
simulations showed that the proposed multi-telescope concept of PLATO can fulfil the defined scientific goal of measuring 
more than 20000 cool dwarfs brighter than my =11 with a precision better than 27 ppm/h which is essential for the study 
of earth-like exo-planetary systems using the transit method. 
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1 Introduction 

The main objective of the PLATO mission is to study exo- 
planetary systems around bright stars using ultra-high preci- 
sion photometry to facilitate ground-based follow-up stud- 
ies. Asteroseismology of the host-stars will be the key to 
pinpoint characteristics such as the age of a system or the 
mass and density of the detected planets. A detailed descrip- 
tion of the PLATO mission can be found in Catala (2009) 
and Claudi (2010). 

The motivation to go to space to acquire photometry of 
stars for asteroseismology or detecting exo-planets is mainly 
the lack of atmospheric disturbances and the diurnal cycle. 
PLATO will be situated in a Lissajous orbit around the La- 
grange point L2 which will lead to a minimization of exter- 
nal influences from Earth. A single field in the sky can be 
monitored for months to years with a very high duty cycle 
and with the same instrument, which leads to a very homo- 
geneous data set, much lower photometric noise levels, and 
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the avoidance of daily aliasing in power spectra. Neverthe- 
less, still many noise sources remain and must be quantified 
in detailed studies to estimate their impact on the quality of 
the observations. The main unavoidable natural noise source 
is photon shot noise, whose relative effect on data can be 
minimized by collecting more photons. The instrumental 
noise sources which are due to the technical construction of 
the spacecraft and due to the limits of the electronics must 
be minimized to ensure that photon noise remains the ma- 
jor contributor of the overall noise budget. Due to the com- 
plex interplay of the noise sources, numerical simulations 
are ideally suited to study the overall performance of space- 
based observations. 

PLATOSim is a software tool dedicated to the realis- 
tic simulation of CCD photometry of stars in the optical 
range acquired by a space-based instrument. It is designed 
in a way to account for the most important natural and in- 
strumental noise sources that have an impact on the perfor- 
mance of the instrument. Due to the need of such a simulator 
for the assessment study for PLATO Mission, we created 
a homogeneous software tool based on pre-existing codes 
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Fig. 1 Screenshot of PLATOSim. On the left panels, pa- 
rameters for the simulations can be set. The right panel 
shows a simulated image which contains several visible 
noise features such as saturation smearing (bright vertical 
spikes of stars) and frame-transfer trails (vertical lines). 



that had been developed for the Eddington and MONS space 
missions (Arentoft et al. 2004; De Ridder et al. 2006). 
PLATOSim has been coded in C++ and features a graphical 
user interface (see Fig. [T]) based on Qt4. The contributions 
of the different noise sources are described through mathe- 
matical models and use a large number of parametrised in- 
puts and numerical pre-computed input data like PSFs, jitter 
time-series, and sky background brightness. 

Several improvements have been applied to the simu- 
lator compared to the pre-existing version. These concern 
mainly a significant increase in computational speed, a more 
realistic treatment of jitter, and the inclusion of two different 
photometric algorithms and statistical tools to estimate the 
noise properties of the data. The graphical interface permits 
in an easy way to test several different instrument set-ups 
by modifying the parameters for the simulations. 

2 Modelling of CCD observations 

The modelling of the synthetic CCD image is carried out 
in several steps which are described in the following. The 
first step is to define the field-of-view (FoV) of one CCD. 
Important parameters here are the size of the CCD in terms 
of the number of pixels, the pixel-scale on the focal plane, 
and the positions and magnitudes of the stars. PLATO will 
consist of a large number of small telescopes, each with four 
CCDs in the focal plane. As a simplification, we only sim- 
ulate the observations acquired by a single CCD and ex- 
trapolate its properties on the complete system by assuming 
that all telescopes behave in a similar way. Moreover, due to 
computational limitations, not the complete CCD but only 
a sub-image of typically 200x200 pixels is simulated. On 
such an image, all stars are assumed to have the same effec- 
tive temperature and the same point spread function (PSF). 
Different distances from the optical axis can be simulated by 
using corresponding PSFs that were pre-computed for dif- 
ferent distances. The positions of the stars on the CCD are 




Fig. 2 Simulated images at sub-pixel scale of a stellar 
field using PSFs at different distances from the optical axis 
(left panel: 0°, middle panel: 6.5°, right panel: 14.1°). Note 
the grid-like structure at the image background which arises 
from the lower intra-pixel sensitivity. Image generated with 
PLATOSim. 

computed assuming orthogonal projection which is a good 
approximation for a relatively small FoV. 

We call the pixels that represent the real pixels of the 
CCD normal pixels. In a sub-image, each normal pixel is 
represented by typically up to 128x128 sub-pixels. These 
sub-pixels assure a realistic modelling of intra-pixel sensi- 
tivity variations of the CCD and of the ACS jitter movement 
which is equivalent to a few percent of a normal pixel. Only 
at the very end of the CCD modelling, the sub-pixels are 
re-binned to the normal pixel scale. After this first step of 
modelling, a noise-free image is created, containing the lo- 
cations and the fluxes of all stars in the field-of-view. 

The next step of the image simulation takes into account 
the physical noise sources of the CCD, the jitter movements 
of the spacecraft, and the instrumental optics. We model the 
global sensitivity variations (flat field) of the CCD by as- 
suming a 1/f spatial power spectrum which resembles that of 
a typical CCD (De Ridder et al. 2006). Additionally, at sub- 
pixel level sensitivity variations that have white-noise char- 
acteristics are introduced. The lower sensitivity between sin- 
gle pixels is also modelled at sub-pixel level. 

Due to various external forces such as the solar wind 
and the magnetic field, and internal disturbances such as the 
reaction wheel noise and structural and thermal flexibility, 
the pointing of the satellite undergoes steady perturbations. 
This jitter has to be kept as small as possible due to its dete- 
riorating effect on the observations by increasing the noise 
of the measured light-curves. We used a jitter time- series 
that has been modelled specifically for PLATO and vali- 
dated using data from two space crafts, ISO and CoRoT, 
where pointing error history is available. The CCD image is 
created by summing up sub-exposures with their individual 
shift of the FoV, which is in the order of a few sub-pixels. 

The image is further degraded by the telescope optics 
which are characterised by the PSF. We used mono-chromatic 
PSFs that have been pre-computed for three different opti- 
cal designs and different distances to the optical axis (see 
Fig. [21) for a wavelength range between 400 and 1000 nm. 
White light PSFs were integrated over the effective wave- 
length range of PLATO taking into account its wavelength 
dependent transmission and quantum efficiency. The result- 
ing image is constructed by convolution in Fourier space. 
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The sky background of each sub-image is assumed to be 
uniformly lit by the zodiacal light and background objects 
like unresolved stars and galaxies. Through photon noise, 
the sky background increases the noise of the stellar photo- 
metric measurements. Values for the sky background are in- 
terpolated from tabulated data from Wehrli (1985) and Lein- 
ertet al. (1998). 

The PLATO CCDs will be working in frame-transfer 
mode to read out the image. This read-out mode causes 
the flux of all stars to be spread along the columns of the 
CCD and introduces an additional noise source by increas- 
ing the background flux. Although we are only simulating 
sub-images of the CCD with PLATOSim, for the determi- 
nation of this structure the stars on the complete CCD are 
taken into account. As it will be done for the PLATO CCDs, 
we save this trailing structure and use it to correct for during 
data reduction. 

In the further steps, the simulated sub-image is re-binned 
to normal pixel scale and the photon noise is computed for 
each pixel. For pixels, where the flux exceeds the full well 
capacity (^saturation), the flux is smeared along the readout 
direction. Finally, charge-transfer efficiency degradation, read- 
out noise, and digital saturation are computed and the image 
is written as a FITS-file to the disk. The series of FITS files 
can be further analysed and reduced with standard scientific 
software. The right panel of Fig. [T] shows the example of an 
image that has been simulated with PLATOSim. 

We also implemented two photometric algorithms in 
PLATOSim, simple aperture and weighted mask photome- 
try, to acquire photometric measurements of all stars on the 
simulated CCD in an efficient way and to compute statis- 
tics to determine their noise levels. Both algorithms take 
into account the sub-pixel position of the stars and the the- 
oretical shape of the PSF to define the photometry mask 
which includes between 60 and 95% of the stellar flux. We 
were able to demonstrate that aperture photometry produces 
lower noise for brighter targets where pollution due to crowd- 
ing in dense fields does not pose a problem. In contrast, 
weighted mask photometry performs generally better for 
fainter stars in a crowded field due to its smaller size of the 
mask. 

3 Simulation study for the PLATO 
assessment phase 

We made an extensive set of simulations to test the perfor- 
mance of the photometric observations of PLATO for dif- 
ferent set-ups of the telescope optics and the CCD in dif- 
ferent positions in the sky. For these simulations, we used a 
catalogue of the proposed FoV of PLATO containing stars 
with my < 15 (M. Barbieri, priv. comm.). Three differ- 
ent sub-fields with spatial stellar densities that account for 
each about a third of the FoV were studied in detail (named 
sparse, medium, and dense field). We computed time-series 
of typically 200 images and applied different photometric 
algorithms for all stars on the CCD and assuming that the 



stars have constant intrinsic magnitude. In a few simula- 
tions, low-amplitude stellar variability has been introduced 
to test the effect of crowding on the detectability of pulsa- 
tions. 

The following questions have been tackled during the 
assessment study of PLATO: 

- How many stars are affected by crowding of other sources 
due to confusion during photometry dependent on the 
number of stars per deg^. Due to the large pixel size of 
15 arcsec/pixel, confusion, i.e. the presence of the target 
star plus at least one fainter (or brighter) star in the pho- 
tometry mask, is a serious problem for stars fainter than 
my = 11. 

- To which extent does confusion influence the detection 
of variable stars? To explore this question we calculated 
how much the photon noise of a target may be elevated 
due to confusion before a pulsation amplitude of 10 ppm 
fails to be statistically significant detectable during a 3- 
year run. Different scenarios, where half or all of the 
polluting sources pulsate with the same amplitude were 
explored. 

- How do different photometric algorithms - simple aper- 
ture and weighted mask photometry - perform in terms 
of the noise budget? Each time- series of simulated frames 
has been analysed with both photometric algorithms. 

- What is the effect of ACS jitter on the noise budget? 
Simulations with and without considering jitter with an 
rms of 0.2 arcsec have been compared. 

- Which optical design performs best? The PSFs of three 
different proposed optical designs and different distances 
to the optical axis were compared. 

- How many stars are observable for PLATO at specific 
noise levels? This is essential to decide if the selected 
baseline design of PLATO is compliant with the mission 
requirements. 



3.1 Results 

Our findings can be summarized as follows. 

- About 80% of all stars with my < 12 in the sparse 
and medium field, which are representative for the major 
part of the PLATO FoV, are not affected significantly by 
crowding from close-by sources. 

- Due to the smaller size of the masks, weighted mask 
photometry is less sensitive due to crowding than aper- 
ture photometry. 

- Stars that are affected by crowding from close-by sources 
also show a larger increase in noise due to ACS jitter. 

- Aperture photometry generally delivers lower noise for 
my < 11, whereas weighted mask photometry per- 
forms better for my > 12 for simulations with and 
without jitter. 

- All three proposed optical designs perform similar with 
respect to crowding and the noise budget. 
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Fig. 3 Expected noise level in ppm/h for observations 
with 40 telescopes using realistic noise modelling without 
jitter in the sub-field and a PSF at the optical axis. Each red 
cross represents the noise of a star using aperture photome- 
try. The dark blue line is the median of the noise in bins of 
0.5 mag. The theoretical photon noise is indicated as green 
line. Measurements that fall below this value are affected by 
smearing from saturated stars. The two relevant noise levels 
for PLATO at 27 and 80 ppm/h are marked. It can be seen 
that the limit for a precision of 27 ppm/h is at approximately 
my = 11. Towards faint magnitudes, the noise increases 
rapidly due to the confusion issue. 



- Applying photometry by using a mask derived from a 
slightly inaccurately modelled PSE delivers no signifi- 
cant signal-to-noise degradation if the used PSE deviates 
by a few degrees from the optical axis from the intrinsic 
PSE or if the stellar colour is not accurately known. 

- ACS jitter significantly increases the measured noise of 
the stars at all magnitudes. The increase is more pro- 
nounced in the dense field and for the near-edge PSEs. A 
subsequent jitter correction for photometry will be nec- 
essary to reach the required specifications of PLATO. 
Such an approach has been proven successful for CoRoT 
(De Oliveira Eialho et al. 2007, Drummond et al. 2006). 

- The number of detectable pulsating stars, assuming pul- 
sation amplitudes of 10 ppm in a broad frequency range, 
decreases when more neighbouring stars add noise to 
the target star due to confusion and/or variability. Still, 
confusion is very limited except for stars in the densest 
fields and for my > 12. 

- The important magnitude range for asteroseismology be- 
tween my = 4 to 1 1 is only weakly affected by confu- 
sion. 

- Star counts reveal that with 40 refractive telescopes (most) 
scientific requirements can be fulfilled. 

- 21000 (goal 20000) dwarfs/sub-giants at 27 ppm/h 

- 13500 (goal 1000) dwarfs/sub-giants with my < 8 
at 27 ppm/h 

- 3 16000 (goal 200000) dwarfs/sub-giants with my > 
8 at 80 ppm/h 




Magnitude 

Fig. 4 Results of simulations to test the detectability of 
a stellar oscillation signal with an amplitude of 10 ppm 
through confusion with neighbouring (pulsating) stars. The 
different coloured bars show the situation for the follow- 
ing assumptions (the target star is assumed to pulsate with 
an amplitude of 10 ppm): Red bars: Only the target star is 
pulsating with an amplitude of 10 ppm. The other stars in 
the photometry mask only contribute noise through photon 
noise. Green bars: 50% of the polluting stellar sources in 
the photometry mask pulsate with an amplitude of 10 ppm 
in a broad frequency range. Blue bars: All polluting stellar 
sources in the photometry mask pulsate with an amplitude 
of 10 ppm in a broad frequency range. Above my = 12, 
no pulsations with an amplitude of 10 ppm can be detected. 
These simulations were done for a medium dense field and 
a PSE at 6° from the optical axis. 



- 52000 (goal 80000) dwarfs/sub-giants at 54 ppm/h 
and later 27 ppm/h 



4 Future outlook 

The ESA Science Program Committee has decided in its 
meeting on Eeb. 17/18, 2010 that the PLATO mission pro- 
ceeds to a Phase B definition study which will end in June, 
2011. A modification of the current design of the mission, 
to enable compliance with the mission requirements regard- 
ing especially the mass budget and schedule is necessary. 
Therefore, in the new design for PLATO, the number of tele- 
scopes will be reduced, the telescope optics will be modified 
(fewer lenses), the EoV will be increased and the observa- 
tion of brighter stars is anticipated. These modifications will 
require more simulations to fine-tune the mission design and 
ensure that the scientific requirements will be met. 
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